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THE  LINEAR  DEPENDENCY  STRUCTURE  OF 
COVARIANCE  NONSTATIONARY  TIME  SERIES  * 
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ABSTRACT:  The  linear  dependence,  feedback  and  causality  structure  of  covariance 
nonstationary  time  series  is  developed.  At  every  instant  in  time,  the  amount  of  linear 
dependence  between  time  series  vectors  is  expressible  as  the  sum  of  the  amount  of  feed¬ 
back  from  the  first  time  series  vector  to  the  second,  the  amount  of  feedback  from  the 
second  time  series  to  the  first  and  the  amount  of  instantaneous  feedback.  The 
parametric  modeling  of  multivariate  covariance  nonstationary  time  series  and  the  com¬ 
putation  of  their  interdependency  structure  from  the  fitted  model  are  also  treated.  The 
time  series  is  modeled  by  a  multivariate  time  varying  autoregressive  (MVTVAR)  model. 
The  fitted  MVTVAR  model  yields  an  instantaneous  power  spectral  density  (IPSD1 
matrix.  The  IPSD  is  used  in  computing  the  linear  dependency  structure  of  nonstationary 
time  series.  An  example  of  the  modeling  and  the  determination  of  instantaneous  causal¬ 
ity  from  a  human  implanted  electrode  seizure  event  EEC  is  shown 

KEYWORDS:  Information  theory,  rime  series,  time-varying  model,  autoregression, 
feedback,  causality,  electroencephalogram. 
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1.  INTRODUCTION 


The  linear  dependence,  feedback  and  causality  structure  of  covariance  nonstationary  lime 
series  is  developed.  The  parametric  modeling  of  multivariate  covariance  nonstationary  time  series 
and  the  computation  of  their  interdependency  structure  from  the  fitted  model  are  also  treated. 
The  analysis  is  applied  to  a  spontaneous  human  epileptic  seizure  event  electroencephalogram 
(EEG).  ' 

The  structure  of  multivariate  stationary  time  series  has  been  of  interest  in  econometrics, 
oceanography,  meteorology,  engineering  and  statistics.  System  quantities  such  as  the  transfer 
function  between  time  series  and  quantitative  measures  such  as  the  amount  of  information  betewen 
time  series  and  the  amount  of  feedback  are  of  interest.  In  the  analysis  of  econometric  data  and 
human  epileptic  seizure  data,  the  presence  or  absence  of  causality  are  additionally  of  interest.  The 
modeling  of  multivariate  stationary  time  series  has  been  treated  for  example  in  Gelfand  and 
Yaglom  1959,  Caines  and  Chan  1975,  Gustafson  ,  Ljung  and  Soderstrom  1977,  Gevers  and  Wertz 
1982  and  Geweke  (1982). 

A  statistical  interpretation  of  causality  in  stationary  time  series  from  nonrepeated  experi¬ 
ments  appears  in  Wiener  1956,  Granger  1963,  Caines  and  Chan  1975,  Geweke  1982  and  references 
therein.  In  stationary  time  series  econometric  data,  driving  or  causality  is  usually  defined  to  mean 
the  presence  of  a  "feedback  free"  condition.  In  stationary  time  series  epileptic  event  EEG  data, 
causality  has  been  defined  operationally  as  the  detection  of  a  delay  of  one  time  series  with  respect 
to  the  causal  time  series,  Brazier  1972,  Gotman  1983.  In  covariance  nonstationary  EEG  data,  the 
objective  of  analysis  is  ihe  estimation  of  the  amount  of  information,  i  he  amount  of  feedback  and 
the  detection  of  causality,  all  at  an  instant  in  lime.  The  inference  of  linear  dependency,  feedback 
and  causality  in  nonstationary  time  series  are  new  topics 

The  determination  of  the  interdependencies  between  the  electrical  activity  at  lifferent  brain 
sites  is  a  subject,  of  current  interest,  in  EEG  analysis  The  recent  work  of  (levins  et  al  1983.  in 
event  related  potentials  (ERPs)  and  of  Mars  and  Lopes  da  Silva  1983  Salto  and  Ilarashima  198J. 
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Innouye  et  al  1983  and  Gotman  1983  in  ongoing  or  background  EEGs  is  evidence  of  that  interest. 

The  literature  on  multivariate  nonstationary  time  series  modeling  is  surprisingly  sparse  when 
we  consider  the  fact  that  many  natural  and  social  science  phenomena  are  nonstationary.  Rosen¬ 
berg  1973,  Sarris  1973  and  Swamy  and  Tinsley  1980  are  studies  in  the  econometrics  literal ure  on 
modeling  time  series  with  time-dependent  stochastic  parameters.  Except  for  Bohlin  1976.  (scaiar 
nonstationary  time  series),  these  methods  have  not  been  applied  to  the  analysis  of  nonstationar> 
EEGs.  The  approach  described  here,  a  deterministic  regression  modeling  of  multichannel  noiista- 
tionary  time  series,  differs  from  the  stochastic  regression  coefficient  modeling  methods  described  in 
the  literature  cited  above. 

The  complex  demodulation  analysis  by  Walter  1968  was  very  likely  the  first  nonstat  lonary 
time  series  analysis  of  the  EEG.  Subsequent  nonstationary  analyses  were  by  Kawabata  1973 
(overlapping  periodograms)  and  Bohlin  1976  and  Isaksson  1975  (nonstationary  time  series  models 
and  Kalman  filter  algorithms).  An  indirect  approach  to  the  analysis  of  nonstationary  covariance 
time  series  is  to  segment  them  into  locally  stationary  time  series  segments  and  to  model  the  seg¬ 
ments  separately  as  stationary  time  series,  Bodenstein  and  Praetorious  1977.  Sagan  and  Sanderson 
1980  and  Benveniste  and  Basseville  1984  and  references  therein.  None  of  the  aforementioned  ana¬ 
lyses  treat  multichannel  EEGs. 

In  Section  2  the  linear  dependency  structure  of  stationary  lime  series,  l  he  linear  dependency 
structure  of  covariance  nonstationarv  time  series  and  the  concept  of  causality  are  treated.  In  Sec¬ 
tion  3  the  modeling  of  covariance  nonstationary  time  series  by  the  MVTVAR  method  is  described 
\n  application  of  MV  TV  AR  modeling  of  human  implanted  eiectrode  epileptic  seizure  event  data 
and  the  detection  of  causality  is  demonstrated  m  Section  4.  Section  5  is  a  Summary  and  Discus- 


2.  LINEAR  DEPENDENCY  STRUCTURE  IN  MULIVARIATE  TIME  SERIES 

2  1  BACKGROUND.  STRUCTURE  IN  STATIONARY  TIME  SERIES. 

The  quantitative  measure  of  the  amount  of  information  between  discrete  random  variables, 
Shannon  1948.  is  the  intellectual  precedent  for  the  quantitative  measure  of  the  amount  of  informa¬ 
tion  in  one  time  series  about  time  series.  The  measure  of  the  the  amount  of  information  between 
two  vectors  of  continuously  distributed  random  variables  X  and  Y.  Shannon’s  amount  of  mutual 
information,  is  given  by, 

-  30  -  30 

c.  I!  fX  y(s.y)ln(fx  y(x.y)  fx(x \fY(y))dxdy.  (2.1.1) 

-  30  -  30 

In  (2.1.1),  fx  y[X.y)  and  f  X{*),f  y{y J  are  respectively  the  joint  probability  distribution  func¬ 
tion  of  the  vector  random  variables  „Y  and  Y  and  the  marginal  probability  density  functions  of  .Y 
and  )  .  Iz  is  a  measure  of  the  amount  of  dependence  between  the  random  vectors  A'  and  1  . 
Equation  (2.1.1)  is  the  negative  of  the  entropy  of  the  "true"  distribution  f\y  with  respect  to  an 
assumed  distribution  fXfy  Also  IT  y  is  the  dissimilarity  or  information  divergence  between  the 
alternative  joint  probability  distribution  f  \  y  and  the  factored  or  independent  probability 
descriptions  fx  p  -  f\Jy  of  the  vector  random  variables  A’.}  Kullback  1 958 .  Formally. 
(2.1.1)  is  equal  to  the  amount  of  information,  on  the  average,  per  observation,  to  reject  the  null 
hypothesis  that  the  random  vectors  X .  Y  are  independent. 

Gelfand  and  Yaglom  1959  computed  (2.1.1)  for  'he  case  of  ,Y.  J  stationary,  jointly  normally 
distributed  •true  series  Let  {;(/)}  be  a  vector  of  time  series  that  is  partitioned  into  the  two  com¬ 
ponent  vector  time  series,  z(t)  [/(().  (/((  I]7".  Let  the  power  spectra!  density  matrix  of  {;(()} 
be  written  in  the  partitioned  form 

*„(/)  >«(/) 

(2.12) 

•V/)  -V/)  ' 
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Alternate  forms  of  the  Shannon-Gelfand-Yaglom  (S-G-Y)  measure  of  the  amount  of  information 
between  the  two  vectors  of  time  series  are,  (Gelfand-Yaglom  1959), 

1/2 

=  f  Ml  S„(/)I  /I  s„(/)l  I  S„(/)l  )df  (2  13) 

'1/2 

1/2 

I.,,  -  f  Ml  S„(/)l  /I  S„(/)-S„(/)S„,(/)s«»(/)l  )<*/ 

-1/2 

1/2 

I,.,  “  -  f  Ml-^t/M 

-1/2 

In  (2.1.3).  .4!  denotes  the  determinant  of  the  matrix  .4  and  the  symbol  '  denotes  the  complex 

conjugate  transpose.  The  first  formula  in  (2.1.3)  mirrors  the  joint  distribution  versus  the  product 
of  independent  distributions  in  (2.1.1).  In  the  second  formula  in  (2.1.3),  the  term  in  the  denomi¬ 
nator  is  the  residual  spectral  density  of  the  time  series  {7(f)}  after  the  removal  by  regression  of 
the  influence  of  (t/(f)}.  If  the  two  time  series  are  independently  distributed,  {y(()}  does  not 
influence  (x(/)},  the  denominator  of  the  second  term  is  identical  to  the  numerator  and  Iz  ,</  is 
zero,  as  expected.  The  third  formula  in  (2.1.1)  is  a  special  case  of  the  amount  of  information 
between  two  scalar  time  series  {7(f)}  and  {y(t)}.  expressed  in  terms  of  the  spectral  coherence  at 
frequency  /.  IV"V(/).  The  last,  formula  In  (2.1.3)  is  particularly  useful  in  "iir  analysis  of  causal¬ 
ity. 

It  should  be  noted  that  the  information  theoretic  I;  and  /y  y  formulas  hold  only  for  the 
case  of  jointly  normal  -ime  series  or  jointly  normal  random  variables.  Otherwise  /  (  Iy  y  ) 

has  an  interpretation  as  a  measure  of  linear  dependence  between  time  series  (  random  variables  ). 

Geweke  1 9* 2  obtains  other  useful  properties  of  the  linear  dependency  structure  between 
jointly  stationary  time  series  To  illustrate  those  properties  we  need  consider  several  alternative 
linear  projections  of  the  tune  series  (7(f))  and  (y(M)  First  we  ssume  that  the  time  scries 
{;(/))  has  an  AH  representation. 
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oo  _  C 

zt  =  zt  -  2  Amzt-m-  E{ zt\  =  0.  «ar[5t]  =  S  =  r  ^ 

m  =  l  ‘"'j/2 


The  terms  Sj2  and  8X6  defined  below.  Define  four  different  linear  projections  of  on  the 

past  of  Xt,  upon  the  pasts  of  Xt  and  yt  upon  the  past  of  Xt  and  the  present  and  past  of  y(  and 
upon  the  past  of  xt  and  the  past  and  future  of  yt.  That  is,  let 


*l.f  ~  xt  2  x2,t  ~  xt  2  ^2,mxt-m  ~  2  E2 .mift-m' 

m  -  t  m  =  1  m  -  I 


l,!  xt  2  ‘^3,nT^t  m  2  .mitt  m* 


x4.t  xt  2  -  m  2  ^i  rr.y> 


ci,ti  z ,  t  • 


Also,  identify  linear  projections  of  {  t/  ( ^ ) }  in  a  similar  manner. 


y \ .t  xt  2  y^.t  '  xt  2  ^  z.^yt  m  2  ^  i.rnyt-w 


y  i,t  xt  2  3,m  yt  -  m  2  ^*3,  miff  m  ■  1/  4 .  t  xt  2  ^  4,mlff-m  2  ^  4  .miff 


'’ar[y..(l  -  -v 


Now.  in  a  natural  way.  define:  (i)  The  amount  of  feedback  from  the  time  series  {  y(  }  to  the 
time  senes.  {  Xt  }  /  ,z,  (ii)  The  amount  of  feedback  from  the  time  series  { Xt  }  to  the  time  series. 
{!/(}■  I  j  ■  ry  and  (>ti)  The  amount  of  instantaneous  feedback  between  the  time  series  {  Xt  }  and  {  y,  } . 


h> 


(2.1.7) 


K  - 


K'. 


t 


r-. 


=  ln(|t„|/|El!|),  /,_  =  lnflfi 


y  \  I  /  I  “'  y2 1 


r  =  ln(|  VI  f\  VI 


From  (2.1.5).  these  feedback  quantities  also  have  specific  interpretations  as  the  amount  of  informa¬ 
tion  in  the  jointly  Gaussian  distributed  time  series  case  or  as  amounts  of  linear  dependencies  in  the 
general  case.  More  emphatically,  in  the  jointly  normal  time  series  case,  feedback  is  information! 

A  theorem  proved  in  Geweke  1982,  (the  measure  of  linear  dependence  between  stationary 
time  series  is  the  sum  of  the  measure  of  linear  feedback  from  the  first  series  to  the  second,  linear 
feedback  from  the  second  to  the  first  and  instantaneous  feedback),  is, 


Theorem : 


^z,y  "  lj  -y  ’  ly  -x  '  ^zy 


(2.1  8) 


with 


II 

s>» 

1 VI  v 

Ml  VI  ' 

1  VD 

»)  !y  -i 

M|t„l 

1  VI  ) 

Ml  V 

IV) 

Hi)  I, 

M|t„l 

i  vn 

M  1  VI 

,  1 VI ) 

"')  Ixy 

i  vi ) 

I»(\t9  *! 

!  V> 

(2  1.9) 


That  is.  each  of  the  amounts  of  linear  dependence  and  amounts  of  feedback  can  be  computed  from 
t  he  residual  matrices  of  i  he  different  project  ions  it  t  fu-  line  series  j-|  ‘  )  and  ;/  I  I  ' 

Proof:  (Similar  to  Geweke  1982  )  f  rom  (2  I  3)-(2  I  7)  we  observe  lhai 

m iviivi  r-;i.  -'d 


i  .* 

* 

Ml  VI  /H(|t|  'VI 

Then. 

Ml V!  i  V  i 

.  * 

✓ 

M'v:  -V.'  i  m  £; 

V'  1 

l  -  V  1 

—  )  b\  mm*"  r\ 

»'  ■ 

f  ■ 

• 

(  n  )  Its  coii-l  r ii <- 1  loti  of  (2.1  *>) 

V  V' 

,  fv  }rT  so 

-  V  - 

** 

to  (  1  1  ) 

*  ’ 

r 

r  * 

• 

r. 

F;, 

‘  jy-Q-'  s-'.s. 

V  _"k  _  -* 

w  hii'li  lead' 


(iii)  follows  by  symmetry  from  (ii). 


(iv)  Follows  from  |  Sl3|  |  “  I  “'I  an<^  s>'mmetry. 


Another  important  result  in  Geweke  1982  is  that  the  feedback  quantities  can  be  additively 
decomposed  in  frequency.  For  example. 


*/  * 


(2.1.10) 


Similar  additive  relationships  hold  for  dependence  and  the  other  feedbacks.  Equality  holds  in 
(2.1.13)  under  mild  technical  conditions  on  the  regression  coefficients  that  are  straightforward  to 
verify,  see  Geweke  1982  for  details.  (The  invertibility  of  the  time  series  model  is  Lhe  issue 
involved.). 

Consider  the  first  and  fourth  projections  of  the  time  series  (j(f)}  in  (2.1.5).  Let  Sx[f),  be 
the  power  spectral  density  computed  from  the  regression  of  (j:(<)}  on  its  own  past.  Let 
6’z|  2  y(/),  be  the  power  spectral  density  computed  from  the  regression  of  (x(<)}  on  its  own  past 
and  the  past  and  future  of  {y(0}  Then,  Iy  ,*(/).  the  feedback  from  {y(f)}  to  (x(f)}  is, 


/„_,(/)  =  Ml  s’,(/)l/|s,( ,.,(/)!). 


(2  111] 


For  example,  when  (y(/)}  does  not  influence  or  feedback  on  |X(M).  (he  spectral  density 
S !  2  (/)  will  be  identical  to  Sx(f).  the  spectral  density  computed  from  { f }  alone.  Then  the 

ratio  in  (2111)  is  one  and  hence  1,^  „ T  (/ )  is  zero,  again  as  expected. 

There  ire  l  wo  important  extensions  of  the  results  m  12  1  8l-i  2.1.  Ill  I)  The  results  hold 
under  conditioning  or  partial  regression  of  the  series  {x(f  ,.(/(/)}  on  another  vector  of  lime  series 
M0J-  Geweke  1984.  2)  TIip  results  hold  for  nonstationarv  time  series  Some  of  those  results  for 
nonstationary  time  --eries  are  discussed  in  Section  2.2.  A  MVTVAH  coefficient  model  which  per¬ 
mits  estimation  of  these  results  is  discussed  in  Section  3. 


2.2  THE  LINEAR  STRUCTURE  OF  COVARIANCE  NONSTATIONARY  TIME  SERIES 


Here  we  extend  the  notation  and  results  on  the  structure  of  stationary  time  series  in  Section 
2.1  to  nonstationary  covariance  time  series.  Assume  an  MVTVAR  representation  for  the  time 
series  z(t),  again  with  z(t)  =  [j:(f)  t/(f)] 


=  Yj  A m,lzt-m  ~  et ’  E\et\  ~  0’  E[ et~k€t]  ~ 


t+k.t 


(2.2.1) 


m  =  I 


with 


V  _ 
—‘t  ~ 


c, 

cl  T,. 


where  St  /  and  Ty 2  t  are  defined  below. 

Define  the  polynomial  operator, 


00 

A{f,t)  =  I  -  Y 1  Am.texP\-~Xmf}-  <2-2-2) 

m  =  1 

Then  in  a  natural  way,  as  an  extension  of  the  definition  of  power  spectral  density  for  stationary 
time  series,  define  the  instantaneous  power  spectral  density.  (IPSD),  matrix  in  terms  of  the 
MVTV  AR  model, 

S(f.t)  -  A(fA)-%A{f,trl*  (2.2.3) 

[The  concept  of  an  instantaneous  power  spectral  density  was  introduced  in  Page  1952.  The  most 
notable  development  since  Page  1952  appears  to  be  by  Priestly  1965.1967.  Priestly  1981,  Chapter 
11  i.»  a  review  of  the  literature.  By  the  way  of  contrast.  Page  and  Priestly  are  both  frequency 
domain  spectral  methods,  while  the  MVTVAR  IPSP  ,s  1  parametric  model  method  ) 

As  we  did  in  the  case  of  stationary  time  series.  identify  linear  projections  of  { ) 


(2.2.4) 


x\,t  ~  Xt  ~  ^  ^  *2  ,t  ~  xt  ^  *^2  m,txt-m  ^  &2 

m=  1  m= 1  m  =  1 

00  30  oc  oc 

x3.t  ~  xt  ~  ^  A  3  mjxt-m  ~  ^  ^  3m,tyt  -  m'  XA,t  ~  xt  ~  ^  ^ 4m.txt-m  ~  ^  -®4  m.tVi-r 

m  =  l  m- 0  m  =  l  m  =  -oc 


™rlx.,tJ  '  -i ,<• 

Also  in  a  similar  manner,  identify  linear  projections  of  {t/;} 


(2.2.5) 


y  1  .r  ::  Vt  S  Xm.tVt  -m~  2/2,1  ~  Vt  ^  ^'2m.txt  -ry  ^  ^2m,.'l/l-m- 

m  =- 1  m  - 1  m  - 1 


oc 

2/3.1  ~  y<  ^  ^ 

m  -  1 


3m, 1  xi  m  ^ 


^  2/4,1  2/f  ^  ^  ^A  m,l!/l-m' 

m  =  0  m-1  m  =  -  0 c 


var(y«,f  ]  = 

Then  using  algebra,  we  obtain  results  for  nonstationary  time  series  that  are  analogous  to  those 
obtained  for  stationary  time  series.  The  amount  of  linear  dependence  at  time  t  between  (j()  and 
(y»).  the  amount  of  feedback  at  time  t  from  (  yt }  to  {£;}.  the  amount  of  feedback  at  time  t  from 
{ X, )  to  {  y(  }.  and  the  instantaneous  amount  of  feedback  at  time  t  from  { J", }  to  {  y( }  satisfy: 

^ i.y.t  ~  ^t  —  y,i  ~  ly  —  z ,1  ~  ^ iy,t  (2.2.6) 
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The  main  result  in  (2.2.6)  is  that  at  each  instant  in  time,  the  amount  of  linear  dependence  between 
time  series  is  the  sum  of  the  amount  of  linear  feedback  from  the  first  series  to  the  second,  linear 
feedback  from  the  second  to  the  first  and  instantaneous  feedback.  Also,  as  in  stationary  time 
series,  in  nonstationary  time  series,  there  is  an  additive  decomposition  of  the  feedbacks  with  fre¬ 
quency.  For  example 

f  Iv-z.t(!)df-  (2-2.7) 

An  example  of  a  relationship  between  the  linear  dependency  measures  for  nonstationary  and 
stationary  time  series  is 

T 

h.y  =  ~  fh.yjt  (2  2.8) 

The  other  nonstationary  linear  dependencies  satisfy  similar  relationships. 

2.3  CAUSALITY  IN  TIME  SERIES 

The  emphasis  of  the  applications  of  structure  in  stationary  time  series  in  Geweke  1982  is  on 
feedback  between  economic  time  series.  As  shown  in  Section  2.2  instantaneous  feedback  (linear 
dependency;  between  nonstationary  time  series,  can  be  expressed  in  simple  formulas.  Indications 
of  ways  of  estimating  those  linear  dependencies  from  nonstationarv  time  series  are  shown  in  Sec¬ 
tion  3.  The  emphasis  here  is  on  a  concept  of  causality  between  time  series.  The  motivation  for 
ihis  section  is  the  determination  of  causality  in  implanted  elect  rode  KKC  data  in  humans  observed 
during  epileptic  seizures.  Our  concept  of  causality  between  time  series  is  different  than  the  "feed¬ 
back  free"  condition  definition  used  in  econometrics  (Geweke  1982  and  references  therein)  and  the 
time  precedence  definition  used  in  tile  analysis  of  epileptic  event  KKG  lata  {Brazier  1972.  Gotman 
19'*.'!)  We  define  causality  in  terms  of  the  S-G-\  and  conditional  S-G-d  measures  between  time 
series  Consider  the  time  series  or  vectors  of  tune  series,  { I  (  t  )  ) .  {(/(/)}  and  {«'(/)}  Assume 
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that  there  is  pairwise  linear  dependence  between  the  series.  That  is, 

Ij  y  ^0,  lz  u,  7^0?  ,  ly  7^0.  (2.3.1) 

Now  assume  that  each  pair  of  time  series  is  conditioned  on  the  excluded  third  time  series  and  the 
following  results  are  obtained. 

4,yi«>  =°-  4.«m»  ~°-  (2.3.2) 

From  (2.3.2),  removing  the  influence  of  (ut(f)}  from  {x(<)}  and  {y(()}  leaves  that  pair  of  time 
series  linearly  related.  Similarly,  removing  the  influence  of  {y|f)}  from  \l{t)}  and  {ut(f)} 
leaves  that  pair  of  time  series  linearly  related.  However,  removal  of  the  influence  of  (J"(l)}  from 
{y( ()}  and  {lt'(f)}  leaves  that  pair  uncorrelated.  Thus  the  series  (x(f)}  uniquely  explains  the 
linear  relationships  between  three  series.  We  say  that  under  the  conditions  in  (2.3.1)  and  (2.3.2) 
that  the  time  series  (x(f)}  is  causal  to  the  time  series  {y(<)}  and  {U)(0}  (  A  mathematical 

model  that  exhibits  the  properties  in  (2.3. l)-(2.3.2)  is  in  Gersch  1972.) 

Epileptic  event  EEGs  tend  to  be  characterized  by  concentrations  of  spectral  energy  in  a  nar¬ 
row  frequency  band  and  causality  or  driving  may  not  be  present  during  an  entire  seizure.  We  are 
thus  motivated  to  exploit  the  instantaneous  additive  decomposition  of  linear  dependence  in  (2.2.6) 
and  adapt  the  definition  of  causality  in  (2.3.1)  and  (2  3.2)  to  a  definition  of  causality  at  some  par¬ 
ticular  frequency  /  j  and  the  time  instant  t . 

Let,  the  IPSI)  matrix  of  the  three  nonstattonary  time  series  (x(/ )  ,y  ( t ) .  tt’  ( t ) }  be  expressed 
:n  component  form. 


s(f.t) 
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1  1 


■  e. 


■W->.dvV-/-V-V-  «•-  z- 
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.1.  .  .IIIIJ  _  !-■  I  M/ 


Several  additional  ingredients  are  needed  for  our  causality  analysis.  The  following  are 


defined:  W ’]y[f-t),  the  instantaneous  spectral  coherence  between  the  two  generic  time  series 

{*(0'V(0)  at  time  t  and  frequency  /,  and  the  instantaneous  partial  spectral  coherence 
S,z  u, (/.<),  Syy  w{f-t)  and  SJy  w(f.t).  respectively  the  instantaneous  spectral  density  of 
conditioned  on  (tr(f)}.  the  instantaneous  spectral  density  of  (t/t/)}  conditioned  on 
{U’(f)},  and  the  the  instantaneous  cross  spectrum  between  { JT  ( f  ) }  and  {y(f)}  at  t,rne  t  and  fre¬ 
quency  /  conditioned  m  (u>(f)J.  In  terms  of  'he  components  of  the  -pertrai  density  matrix  in 
(2.3.3)  these  are  given  by 

ir^/.r  -  (/.*)!  1% i /■•'> 

u.(/.n  -  H;u, (/./)) 

*yy  «(f-t '  -  »•;,(/•/)) 

Sxyiw(f,t)  -  S:y(f.t)Sy:(f. t)  SwJf.t) 

The  instantaneous  partial  spectral  coherence  between  the  generic  time  series  {/■(?). y(?)} 
conditioned  on  the  time  series  {?/'(/)}  at  time  t  and  t  and  frequency  /.  computed  from  the  instan¬ 
taneous  partial  spectra  in  (2.3.1)  is 


J/-0I 


i  s 

■  yy  to 


Now.  using  the  additive  frequency  domain  decomposition  property  of  feedback  between  time 
series  (2.2  7).  and  our  definit  ion  of  causality  in  (2  3. 1 )  and  (2  3.2)  we  define  i  he  senes  {/(/)}  to  be 
causal  in  i  lie  .Ties  ;  t/(  /  ! . //•  l  /  I  n  freouem  >  ]  (  h  lme  !  if  t  fie  following  'onditions  are 
sat  istied 


/..  ,.•(/,  I  -n.  7. ,,  ,(/.,)  7v ,  .(/, )  •  i).  (2  3  ii| 

/•  ,  '  ■  <>-  yt'/ti  "•  /vvr.;(/-ll  "• 

(•  rorri  the  thirl  equal  e  in  in  (2  I  3 )  we  seen  I  fiat  I  fie  defi  n  u  ion  of  causality  in  j  2  3  *>l  is  equ  i\  a  lent 
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to  the  following  conditions  on  spectral  coherences  and  partial  coherences 


r  o.  WljfAJ)  /  0,  WIMa'O  *  0-  (2  3  7) 

Ky  Ma- 0  -  0.  Wlw  y(fA,t)  ,  0.  W'y  w  A/a-0  =  0. 

Thus  the  identification  of  causality  at  frequency  f  A  at  time  t  is  determined  by  the  detection  of 
zero  partial  instantaneous  coherence  This  definition  of  causality  at  an  instant  in  time  is  a  general¬ 
ization  of  a  related  concept  for  stationary  time  series  introduced  in  Gersch  and  Goddard  1970  and 
used  in  Gersch  1972  Gersch  and  Tharp  1976  and  Bnllinger  1976 

The  list ribut ion  of  coherence  and  partial  coherence  is  'reared  in  Hannan  1970.  Flnilinger 
1974  and  Koopmans  1974  For  our  purposes  the  most  important  result  is  the  distribution  of  a 
transformed  version  of  the  partial  spectral  coherence  in  the  vicinity  of  zero  partial  coherence  A 
convenient  form  of  that  result,  Koopmans  1974.  is  given  by 


W-1)K  A/} 


[1  -  K  A/)) 


F ?  21  IS-  1)' 


(2.3  8) 


In  (2.3.10).  U  is  the  number  of  degrees  of  freedom  in  the  F  distribution  with  2  and  2(i/  1) 


degrees  of  freedom  in  the  numerator  and  denominator  respectively 


3.  A  MULTIVARIATE  TIME  VARYING  AR  COEFFICIENT  MODEL  AND 
LINEAR  DEPENDENCY  COMPUTATIONS 

A  parametric  model  is  assumed  for  the  time  series.  The  successful  use  of  autoregressive  (AR) 
models  in  stationary  time  series  analysis  motivates  consideration  of  a  time-varying  AR  coefficient 
model  for  nonstationary  covariance  time  series  modeling.  Denote  the  D  component  row  vector  at 
time  n  by  r(n)  =  ( r, ( n ) , ...,rD ( n) )  T .  Then,  the  multivariate  time-varying  AR.  (MVTVAR). 
model  is 

z(n)  -  4,„z(n-l)  -  a2nz(n-2)  -...-  -4p„x(n-p)  -  <(n).  (3.1  1) 

£(«(«))  =  0,  £(<(n  +  *)t(n)T)  =  V(a)6n.tn 

In  (3.1  1),  the  { .4,  „i  =  l,...p;  n  =  l....,.V}  are  DxD  coefficient  matrices.  The  pxDxDx A’  unk¬ 
nown  AR  model  parameters  and  the  \xD(  D  ~  \  )  2  parameters  in  V'(n).  n  1 . V  in  the  model  in 

(3  11)  are  to  be  estimated.  There  are  more  parameters  than  data,  so  least  squares  or  maximum 
likelihood  methods  for  estimating  the  unknown  parameters  will  not  yield  useful  results  The  unk¬ 
nown  parameters  can  be  modeled  implicitly.  A  strategy  for  economizing  on  the  number  of  param¬ 
eters  to  be  estimated  is  to  consider  the  vector  of  time  series  one  component  at  a  time  and  to  regress 
that  time  series  upon  a  lagged  version  of  itself  and  upon  t  he  other  components  of  the  sector  of 
time  series  in  ari  orthogonal  polynomial  least  squares  method  of  modeling  That  is.  express  each  of 
the  elements  in  the  matrix  of  time  varying  AR  coefficients  as  a  linear  combination  of  say  J  orthog¬ 
onal  polynomial  functions  ,(  time,  where, 1  is  a  small  number  orupare.j  in  \  'tie  number  of  obser¬ 
vations  I  sing  'his  method  as  many  is  p\l)\J  .efficients  are  titled  to  each  of  the  i>  component 
time  series  If  tills  number  is  small  compared  to  A.  a  reasonable  model  can  be  fitted  The  total 
number  of  fitted  coefficients  in  the  model  is  then  pxDxHx.l  I  hat  number  is  considcrablv  smaller 
than  the  number  of  implicit  \R  coefficients  m  1  It  A  The  \|\'TV  \R  model  was 

introduced  m  tiers,  h  and  Kitagawa  TtM.3  More  spec  di<  all\  lei 
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A,  »  =  'a,*(«.n)|.  r,k  =  l....,D 


(3.1.2) 


J  \ 

art(«,n)  =  Y]c,*(t,»/(;,n),  *  =  1,-  ,P 

i-o 


.V 


/(>»  =  S(-i)* 
•  =  0 


n'*>  A’1*1. 


J 


n  -  0. 1 . V 


In  (3.1.2),  art(t,n)  r,k  =  l,...,D  are  the  elements  of  the  time  \arying  AR  coefficient  matrix 
.4,  „  i  =  1  ...p;  n  =  1 N.  Also  in  (3.L.2)  =  n'  '(n  — *)!.  .V^'1*  ■=  .V!  !'  where  n!  is  n  fac¬ 

torial  and  N  is  the  number  of  data  points.  The  functions  f(j.n)  that  we  use  in  (3.1.2),  are  the 
discrete  orthogonal  Legendre  polynomials. 

The  orthogonal  polynomials  satisfy 

,v 

£/(j,n)/(*.n)  ^  0,  k*j.  ('■  ‘ 

n  =  0 

The  first  three  discrete  orthogonal  Legendre  polynomials  are 

/(0,n)  1,  /(l,n)=l-2n/.\\  /(2.n)  =  1  -6n/  .V  -  6n(  n  -  1 )  .V(.Y  1 ) .  (3.1.4) 

From  (3  12.1),  fitting  the  V1VTVAR  model  (3  1.1).  with  the  zeroth  degree  Legendre  polynomial  is 
equivalent  to  fitting  a  stationary  coefficient  multivariate  autoregressive  (MV.\R)  model 

Orthogonal  polynomial-least  squares  modeling  of  scalar  time  series  using  t irne- vary ine  \R 
coefficient  models  appears  in  Kozin  1077  and  Orenier  1083.  (lersch  ami  Kitagawa  1983  is  an 
extension  of  that  method  to  the  fitting  of  multivariate  time  series  using  multivariate  time-varying 
AR  coefficient  (MVTVAR)  models  to  economic  data  in  that  modeling,  the  orthogonal  polynomial 
least  squares  computations  were  realized  ,n  a  Householder  transformation  algorithm.  Xkaike- 
AIC  statistic.  Akaike  1974.  orthogonal  polynomial  model  degree  and  regression  subset  selection 
( r or  more  details  on  llie  modeling  method  see  (lersch  and  Kitagawa  I  '8.3  j 

hach  <t  the  L,  ,  ami  ,  i  I  ,1  quantities  that  are  require. i  ■  ..mpme  in,  n  ,i  m 
laneous  linear  dependency  relations  in  (2  3  ( i )  are  c, input, able  |,\  tiiimi'  ,f  \1\  I  \  \H  m,  .bi- 
Flint  is.  ,  7  j  ,  are  obtained  from  the  fit  ..ft  he  \1\  |\  \((  model  ;.  the  .  ■  f  i  /  ( ]  y  ( f  ] 


j  :> 


Then  ty,  =  t2l  -  CtT2\Cf.  From  (2.2.6),  |  £4  ,|  =  t2M  \  ,  ,|  T,  ,|  The  terms  F,  „  T,  ,  ran 
be  obtained  by  the  fit  of  MVTVAR  models  to  r(f)  and  y(t)  respectively,  t'sing  these  quantities  in 
(2  2  6)  each  of  the  terms  It  ,,  /,  ,.  /#  .  and  /  ,  ran  be  computed. 

Also,  the  1PSD  ran  be  estimated  by  adapting  the  definition  in  (2  2  21-12  2.3)  to  the  finite  lag 
order  fitted  MVTVAR  model  b<>r  the  fitted  \1V'T\  111  model  of  lag  inter  1/,  iefine  the  polyno¬ 
mial  operator 

v 

M f  t)  I  V  I,  , e x p  2r.  tnf 

-n  1 

Then,  define  the  estimated  instantaneous  power  -iie-trai  letisitv  tnatr.x  ,:i  erm 
MVTVAR  model. 

>'(/•')  ii/  n  "L,  u/.ti  1  *  i3i  :,) 

The  spertral  roherenre  and  partial  spectral  coherence  at  frequency  /  at  the  time  instant  I  ran  be 
computed  using  the  estimated  Il’bP 


(X  1  4) 

I  !  fie  fit  I  e<  1 


4.  AN  EXAMPLE,  ANALYSIS  OF  A  HUMAN  EPILEPTIC  EVENT  EEG 


We  illustrate  a  case  in  which  an  epileptic  focus  was  located  and  determined  to  be  present 
during  a  short  time  interval  from  an  analysis  of  an  epileptic  EE(1  event  from  deeply  implanted 
electrodes  in  a  human.  (About  of  individuals  with  epilepsy  have  seizures  initiated  from  an 

anatomically  localized  brain  region  or  seizure  focus.  The  focus  initiated  seizure  propagates  through 
the  brain,  Hauser  and  Kurland  1975  Many  of  these  individuals  do  not  respond  to  drug  treatment. 
In  such  individuals,  if  a  unique  anatomical  site  from  which  the  seizures  emanate  can  be  localized  to 
an  operable  site,  they  may  be  suitable  candidates  for  surgical  treatment  to  remove  the  seizure  focus 
and  potentially  relieve  the  individual  from  seizures.  Anderman  1987  )  Figure  1  illustrates  data 
from  a  7-second  fi  site  data  record  of  the  electrical  activity  at  the  dramatic  onset  of  a  spontaneous 
seizure  event.  (This  data  was  supplied  by  Jeffrey  I*  Lieb.  I  Cl. A  Heed  Neurological  Research 

('enter  The  original  data  w  obtained  at  a  200  sample  second  rale  The  data  used  for  analysis 

was  rate  reduced  to  50  samples  per  second  )  The  objective  of  the  analysis  was  to  determine 
whether  any  one  of  the  observed  anatomical  recording  sites  could  be  interpreted  as  "driving"  or 
being  causal  to  the  electrical  activity  at  the  other  recording  sites  The  data  is  complex  and  there 
is  no  obvious  visual  clue  thai  might  unambiguously  identify  an  .nutating  or  driving  site  <>r  identify 
w  hen  driving  is  present 

I  he  appearance  of  t  he  individual  I.K(I  traces  m  Figure  I  reveal-  that  i  tie  o  a  |  structure  .if 
the  time  series  changes  w.  1 1  h  tune  |'he  analysis  involve-  the  titling  f  a  \1\  l\  \R  model  t . .  i  tie 
data  and  a  subsequent  spectral  analysis  using  the  ll’SI) 

1  1 1  e  :  a  I  a  a  n  1 1  v  sis  will  |>e  .lesvn  la*,  i  hi  t  iriie  oi  i  he  i  n  a  I  v  .  i-  t  ■  n  -  . ,  m  u  it  a  tie*  o-  (  |  (  ,  -  ,  n 

channels  !at  e|e,j  |  2  an)  in  figure  1  In  figure  2  we  -h-  w  a  'liri'-oe"  view  ,.|  i  tie  estimate! 

Il’-I  l  -  versus  frequency  . imputed  at  successive  |  2  '•••  .  nd  o.tervai-  from  i  tie  hue  f  M\  f  \  \  |; 

•no  !e|  .-a  -  !i  -  f  .  Iiatine|s  f  2  and  ?•  Initially  there  are  rel.ative.'v  -h  »r :  u  ,  •  r  .  |  .  . . 

n  -a  (i  .|  .flannels  |  2  and  at  7  25  Hertz  and  It  5  1 1  e  r '  /  'to  1  |  (|  .v  \\ ,•  •  j.,  ,« 

'e  lie  |  It:  '  I  II  idled  HI  atllplll  llde  llld  the  dlstribll'  loti  '  f  speetra.  er  ef|U  lie  ;  i  •  r  il't  I  . 


Later  in  the  record  the  emergence  of  a  spectral  peak  at  4  25  Hertz  can  be  observed  in  each  channel. 
A  current  practice  is  to  analyze  epileptic  event  EEC's  of  no  less  than  6.25  second  duration  in 
overlapping  2.5  second  intervals  by  classical  windowed  periodogram  spectral  analysis  methods  as  if 
the  time  series  in  each  2.5  second  interval  were  stationary,  Gotman  1983, 1987,  Lieb  et  al.  1987. 
(The  windowed  periodogram  methods  used  by  Gotman  and  Lieb  do  not  have  the  the  time- 
frequency  resolution  properties  necessary  to  capture  the  transitory  characteristics  of  the  rapidly 
changing  epileptic  event  KEG  data.) 

We  identify  causality  or  driving  at  an  instant  in  time  by  a  frequency  domain  analysis  using 
the  estimated  1PSD  via  the  concept  described  in  Section  3.  Computational  results  of  both  the 
evolving  spectral  coherence  and  evolving  partial  spectral  coherence  for  each  of  the  three  possible 
pairs  of  three  different  data  channels  at  successive  one  second  time  instants  are  shown  in  Figure  3. 
The  computational  results  are  used  to  determine  whether  and  when  one  channel  drives  the  other 
t  wo. 

A  view  tlown  the  columns  of  the  one  second  apart  "snapshots''  in  Figure  3  illustrates  the 
changing  with  time  structure  of  pairwise  coherence  versus  frequency  between  pairs  of  channels  and 
the  corresponding  changing  with  time  partial  coherence  versus  frequency  At  one  and  two  seconds 
into  the  record,  there  are  sharp  coherence  peaks  between  each  of  the  >f  the  channel  pairs  in  the 
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We  conclude  that  the  pairwise  coherence  between  the  three  channels  is  a  consequence  of  the 


fact  that  channel  1  was  driving  channels  2  and  3  for  the  first  2  seconds  of  the  record  in  the  vicinity 
of  7.25  Hertz.  The  pairwise  coherences  and  partial  coherences  at  three  seconds  and  subsequently 
reveals  a  decrease  in  the  pairwise  coherence  between  each  of  the  channel  pairs  and  cessation  of 
driving  by  channel  1  of  channels  2  and  3.  Later  in  the  record,  the  locations  of  the  peaks  of  the 
spectral  coherence  shifted  to  a  lower  frequency  and  the  partial  coherence  does  not  indicate  driving. 

The  relatively  low  partial  coherence  between  channels  2  and  3  partiaied  on  hannei  1  throughout 
the  epoch  is  compatible  with  evidence  cited  by  Gotman  1987 

f  igure  1  is  a  birds  eye  view  of  the  evolution  of  the  pairwise  coherences  and  the  partial  coher¬ 
ences.  This  illustration  was  computed  from  the  same  data  as  Figure  3  [n  general,  partial  coher¬ 
ences  are  smaller  than  coherences.  The  relatively  flat  partial  coherence  of  channels  2  and  3  oar- 
tialed  on  channel  1  suggests  that,  this  kind  of  illustration  may  be  a  useful  diagnostic  aid  to  identify 
candidate  driving  channels  (Gotitian  1983.1987  is  an  alternative  method  of  detecting  driving 

For  the  purposes  of  comparison,  in  Figure  5  we  show  the  coherence  and  partial  coherence 
results  computed  from  a  multivariate  \R  model,  as  if  the  time  -enes  were  stationary  Such  ana¬ 
lyses  yield  a  blurred  version  <,(  the  sharp  time-frequency  resolution  features  available  from  an 
\1  \  1  \  VK-M’Sl)  analysis  Mso.  these  results  do  not  definitively  implicate  channel  1  as  (driving  the 
a  her  t  wo  hannei- 
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5.  SUMMARY  AND  DISCUSSION 


Measures  of  linear  dependence  and  feedback  at  each  instant  in  time  for  multiple  covariance 
nonstationary  time  series  have  been  proposed  Several  new  ideas  are  introduced.  The  development 
is  heuristic  For  jointly  normally  distributed  stationary  time  series,  linear  dependence  is 
equivalent  to  the  Shannon-*  ieifand- Yaglom  measure  .if  the  amount  of  information  between  time 
series.  The  concept  of  information  at  an  instant  in  time  appears  to  be  new.  The  time  series  are 
assumed  to  be  represented  by  i  multivariate  time  varying  autoregressive  ^MYTVAR)  model.  The 
MYTVAR  model  is  the  key  to  the  computation  of  linear  dependency  and  linear  feedback  at  each 
instant  in  time  At  each  instant  in  time  the  measure  .if  linear  dependence  is  the  sum  of  the  meas¬ 
ure  of  linear  feedback  from  the  :irst  series  to  the  second,  inear  feed  bar  k  from  the  second  series  to 
the  first  an<l  instantaneous  feedback  The  measures  of  linear  feedback  from  one  senes  to  another 
can  be  addittvely  decomposed  by  frequency. 

The  time  evolution  of  AR  parameters  in  the  M VT V  A R  model  is  expressed  as  linear  combina¬ 
tions  of  discrete  Legendre  orthogonal  polynomial  functions  of  time.  The  MVTVAR  model  is  fitted 
by  a  Householder  transformation- Akaike  AIC  method  The  MYTVAR  model  is  exploited  to  intro¬ 
duce  the  concept  of  an  instantaneous  power  spectral  density  (IPS!)).  This  is  done  in  a  natural  way 
as  an  extension  of  the  power  spectral  density  computed  from  multivariate  AR  models  for  station¬ 
ary  time  series  File  •  line-frequency  resolution  properties  of  i  he  M\  T\  AR-IPSI)  computations  are 
-harper  than  those  obtained  by  segmenting  nonsl at tonary  time  series  into  successive  stationary 
segments  and  using  windowed  period. igram  spectral  analysis  methods  .r  by  Priestly's  evolutionary 
peei  ra  met  hi  d. 

I  tie  decomposition  ot  linear  dependence  is  an  analysis  •  I  variance,  tut  a!  tins  stage.  ,a  theory 
of  inference  lias  not  been  developed  Potentially  i  li>-  \!\  1  \  NR  and  I  PS  I )  can  be  useful  for  investi¬ 
gating  'tie  trist  an  t  aneoijslv  changing  mterrei.nl  n  ri  -  h  1 1  >*-  oi  n  ■■  uno.'r  apliic  meteorological  and  other 
multivariate  times  senes  phenomena  as  well  e  invest  ig:n  i.  ms  into  the  nature  of  the  instantane¬ 
ously  changing  ititra  and  inter  reiebral  propagation  of  epilepsy 


A  new  concept  of  causality  at  an  instant  in  time  was  proposed.  A  related  definition  of 
causality  for  stationary  time  series  is  different  than  the  familiar  "feedback  free"  Granger  causality 
in  econometrics  or  the  time  precedence  concept  by  Brazier  in  experimental  neurophysiology.  This 
concept  of  causality  at  an  instant  in  time  appears  to  be  particularly  well  suited  for  the  determina¬ 
tion  of  driving  or  causality  in  epileptic  event  EEGs.  To  date,  the  results  of  epileptic  focus  location 
obtained  by  our  method  have  been  consistent  with  analyses  done  earlier  and  with  the  successful 
outcomes  of  -urgical  temporal  lobectomy  seizure  elimination  procedures.  The  perspective  or 
birds-eye  views  of  coherence  and  partial  coherence  versus  frequency  and  time  appears  to  be  a 
potentially  useful  diagnostic  for  the  identification  of  candidate  driving  sites.  We  anticipate  doing 
computations  on  other  epileptic  event  data  sets  in  humans  and  animals  and  comparing  results 
obtained  by  other  methods. 
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LEGENDS 


FIGF  RE  l.  Intracerebral  EEC  during  a  spontaneous  7  second  seizure  episode.  From  the  top  to 
bottom  the  data  channels  are:  right  (rt.)  rt.  amygdala,  rt.  ant.pes  hippocampi,  rt  post,  pes  hippo¬ 
campi..  rt,  ant  parahippocampal  g>rus,  rt  mid.  parahippocampal  gyrus,  rt.  parahippocampal 
gyrus.  Analysis  of  the  EEG  on  the  left  side  of  the  brain  did  not  indicate  very  substantial  involve¬ 
ment  in  epileptic  activity  in  this  epoch  and  is  not  -hewn  here.  The  liscussion  in  the  'ext  s 
confined  to  the  analysis  of  the  channels  marked  1,2  and  3  in  Figure  1 

FKil  RE  2.  Instantaneous  power  spectral  densities  in  decibels  of  channels  1.2  ana  3  versus  fre¬ 
quency  and  time  in  Hertz  at  successive  1  2  ,econd  intervals 

FUJI  RE  3.  Instantaneous  coherence  (solid  line)  and  instantaneous  partial  coherence  (dotted  line) 
versus  frequency  in  Hertz  at.  successive  1  second  intervals  for  7  seconds  The  number  pairs  before 
each  column  of  graphs  indicate  the  channel  pairs-  The  number  triples  indicate  the  partial  coher¬ 
ence  between  the  first  pair  of  numbers  with  the  influence  of  the  third  channel  removed  from  the 
first  two  channels  by  regression. 

FIGl  RE  4.  Instantaneous  coherence  and  instantaneous  partial  coherence  versus  frequenr;  tr 
Hertz  and  time  at  successive  1  second  intervals  for  7  seconds. 

FIG  I  RE  Coherence  and  partial  coherence  [dotted  lines!  •  ■mputed  is  if  the  time  -cries  acre 
stationary.  The  number  pairs  before  each  column  of  graphs  indicate  the  channel  pairs  The 
number  triples  indicate  the  partial  coherence  between  the  first  pair  >f  numbers  wrh  the  influence 
of  the  third  channel  removed  from  the  first  two  channels  b\  regression 


LEtiENDS 


FIG  IRE  1  Ir:t  racerebral  KEG  during  a  spontaneous  7  second  seizure  episode  From  the  top  to 
bottom  t  he  data  channels  are  right  (rt.)  rt  amygdala.  rt  ant  pes  hippocampi  rt  post  pes  hippo¬ 
campi..  n  mi  paranippocampal  gyrus.  ri  mid  parahippocampal  g\  ru».  rt  paraiuppoc.ampai 
gyrus  Analysis  of  the  KEG  on  the  left  side  of  the  brain  did  not  indicate  ver>  substantial  involve¬ 
ment  ,n  epiietitir  activity  .n  this  epoch  and  s  not  -hown  here  The  iisnission  n  he  eat  s 
confined  to  the  analysts  of  the  channels  marked  1.2  and  3  in  Figure  1 

HGI  HE  2  Instantaneous  power  spectral  densities  in  decibels  of  cnannels  1  2  and  3  versus  fre¬ 
quency  aim  tine  n  Ib'rtz  at  -uc-essive  1  2  second  mtervais 

FIGl’RE  3  Instantaneous  coherence  (solid  linei  and  instantaneous  partial  coherence  (dotted  line) 
\ersus  frequence  .n  Hertz  at  successive  I  second  intervals  f>r  7  seconds  The  number  oairs  before 
each  column  of  graphs  indicate  the  channel  pairs  The  'lumber  triples  indicate  i  he  partial  coher¬ 
ence  between  the  first  pair  of  numbers  with  the  influence  of  the  third  channel  removed  from  t  h  •' 
first  two  channels  by  regression 

HGI  RE  F  Instantaneous  coherence  arid  instantaneous  partial  coherence  versus  rre.;uetic\  n 
Hertz  and  time  it  successive  1  -er  >nd  intervals  for  7  second' 
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